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Either ATIC or Fermi-LAT data can be fitted together with the PAMELA data 
by three components: primary background ~ E^'^'^^ secondary background ~ E~^'^^ 
and an additional source of electrons ~ i?"'''"Exp(— i?/£^cut)- We find that the best 
fits for ATIC + PAMELA and for Fermi + PAMELA are approximately the same, 
7(j ~ 2 and E'cut ~ 500 GeV. However, the ATIC data have a narrow bump between 
300 GeV and 600 GeV which contradicts the smooth Fermi spectrum. An interpre- 
tation of the ATIC bump as well as the featureless Fermi spectrum in terms of dark 
matter models and pulsars is discussed. 
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The question of interpretation of Fermi-LAT [T], HESS [21 |3], and ATIC data [1] can be 
split into two parts: general properties of the flux and the presence of features. If one takes into 
account PAMELA data [S] , then both Fermi and ATIC require the existence of an additional flux 
of electrons and positrons complementary to the standard primary and secondary backgrounds. 

We will consider the following form for the additional flux 

F<j~E"^"e"^. (1) 

The general properties of the flux will be parameterized by the index 7^ and the exponential 
cutoff Ecuf Significant deviations from this form will be considered as "features". 

We will assume the primary background ~ E~'^p, the index 7p ^ 3.3 can be estimated from 
the electron injection index 2 — 2.5 due to shock acceleration in the supernovae explosions [HI E] 
(a similar estimation for the shock acceleration in gamma ray bursts can be found in [HI |9]). 
The shift of the index from 2 — 2.5 to 3.3 is due to cooling during propagation [TH]. Secondary 
background of electrons and positrons is produced by collisions of high energy protons and nuclei 
~ E'^^-^ [Illll2| with the dust. We will assume the secondary background ~ E~"'% 7^ ^ 3.6 where 
the difference in the indices is again due to cooling. Since the normalization and the indices of 
the backgrounds are known only approximately, we will treat them as independent parameters 
in the fits. 

In order to find the properties of the additional flux, we will use the ATIC data, the Fermi 
data, and the PAMELA data above 10 GeV (the points below 10 GeV are assumed to suffer from 
solar modulation and we will discard them for the purposes of current analysis). The results of 
the fits to ATIC + PAMELA and to Fermi + PAMELA are presented in Figs, [l] and |2] and in 
Table |Tj Both ATIC + PAMELA and Fermi + PAMELA are fitted best by an additional flux 
with 7a 2 and Ecut ~ 500 GeV, i.e., these experiments are consistent with each other from the 
point of view of general properties of the flux parameterized by Eq. ([T]). It should be noted that 
without PAMELA data, the ATIC bump is better fitted with a harder additional flux, 7^ ^ 1.7, 
while the Fermi data are consistent with the primary background 7^ ~ 3.0 — 3.1 [H [13] . 

The second question is the presence of features. The Fermi data are well approximated by the 
additional source in Eq. ([T]) while the ATIC data have a deviation from ([T]) between 300 GeV 
and 600 Gev. Thus, although the parameters of the additional flux in Eq. ([T]) necessary to fit 
ATIC and Fermi data are similar, the ATIC data have a bump while the Fermi data don't have 
any features, i.e., there is a disagreement between the two data sets at energies 300 - 600 GeV. 
As one can see in Table [I], the best fit for Fermi + PAMELA has a better reduced chi-squared, 
Xr = 0-4, than the best fit for ATIC + PAMELA, Xr = 1-5- 

In the following we will summarize the properties of dark matter (DM) models and pulsars 
necessary to reproduce both the fiux with features and the fiux without features. Due to energy 
losses, the source of high energy electrons should be close to Earth, within approximately 1-3 
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FIG. 1: The results of fitting the primary and secondary backgrounds together with an additional source 
of electrons and positrons to Fermi + PAMELA and to ATIC + PAMELA. As everywhere else in the 
paper, for PAMELA, only > 10 GeV points are used. The best fit values and the parameter ranges 
are given in Table |l) The star (cross) represents the best fit to Fermi + PAMELA (ATIC + PAMELA) 
with the reduced chi-squared x?best ~ ^-^ (^rbcst ~ The dashed (solid) contours correspond to 

a:? best + 1 and Xr best + 2 for Fermi (ATIC) + PAMELA. The concordance model has = 0.9 (2.0) for 
Fermi (ATIC) + PAMELA. 

kpc. For the purposes of flux calculation, the DM distribution can be viewed as homogeneous 
and constant [H] (unless there is a significant contribution from a local DM substructure such 
as a clump [HI [T51 [IS])- The flux from a homogeneous source is [TU] 

where E = —b{E) is the energy losses. At > 10 GeV the energy losses are due to Inverse 
Compton Scattering and synchrotron radiation in the galactic magnetic field, thus b{E) = BqE"^ 
[TU] . Qbm is the source function for e+e" produced by annihilating or decaying DM, Qbm = 

dN 
dEdVdt ■ 

The cutoff energy ii^cut is the energy where the integral on the right hand side of Eq. (j2]) is 
saturated. For energies E -C -Ecut the integral is insensitive to the variations of the lower limit 
and -Fdm ~ ^/b{E). Thus an index 7^ «i 2 is a universal prediction of DM. Local clumps of DM 
result in a harder spectrum 7^ < 2 [TU [15]. The dependence on the host halo profile leads to 
a softer spectrum 7^ > 2 at low energies p3]- In models with small diffusion height, there may 
be some hardening 7^ < 2 at low energies due to leakage of electrons from the Galaxy (see, e.g., 
model M2 in Figure 3 of [H]). 

For a given E'cut ~ 500 GeV, the DM mass and the shape of the spectrum depend on the DM 
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(GeV) 


(GeV^m-^g-igj,-!) 


X? 


7p 


7s 


ATIC + PAMELA 


1.95 ±0.15 


480 ± 200 


50 ±8 


1.5-2.5 


3.27 


3.63 


Fermi + PAMELA 


2.15 ±0.10 


500 ± 150 


44 ±6 


0.4 - 1.4 


3.31 


3.63 


Concordance model 


2.05 


450 


47 


1.3 


3.29 


3.63 



TABLE L Numerical values for the fits presented in Fig. [T] The normalization is given for Eq = 100 
GeV. The ranges of parameters correspond to the ranges of reduced chi-squared Xr the fifth column. 
The error bars for Fermi are computed as square root of systematic plus statistical errors squared. 7p and 
7s are the indices of the primary and secondary backgrounds respectively. The for the concordance 
model is computed using Fermi, ATIC, and PAMELA (>10 GeV) points. The best fits are found by 
varying 7 parameters: 2 indices and 2 normalization constants for primary and secondary backgrounds 
together with the index, the normalization and the cutoff of the additional flux in Eq. ([T]). 

model. In general, DM annihilation is followed by a sequence of decays leading to electrons and 
positrons together with other stable particles in the end. Models with many steps in the decay 
process have smooth e"'"e~ spectrum and a large DM mass Mdm ^ -E'cut [IH HE]- These models 
are favored [191 120] by featureless Fermi spectrum. Models with decay channels through W , Z 
gauge bosons and quarks have stricter constraints due to absence of significant deviations from 
the expected backgrounds for anti-protons [21] and diffuse gamma rays from the Galactic center 
\Z2\ [23]. The flux of gamma rays from DM clumps assuming hh annihilation channel was also 
estimated [21]. DM models with small DM mass Mdm ~ 500 GeV (e.g., [23 ES]) seem to be in 
tension with Fermi data due to absence of significant step-like features below or around 500 Gev. 

DM models with few decay steps have a sharper cutoff and Mdm ~ -E'cut [H HSl I2Z1- These 
models give reasonable fits to the Fermi data for Mdm ~ 1 TeV [13], however, due to a sharp 
cutoff near Mdm they, generally, fit better the ATIC data [18] . These models may have additional 
constrained due to final state radiation from the Galactic center and the Galactic ridge [iHl [27] . 

Local clumps may produce additional features at high energies. The presence of a large local 
clump is disfavored by Fermi but consistent with ATIC [13] • Furthermore, in models with many 
decay steps, a significant contribution from a local clump may be necessary to fit the ATIC bump. 

Thus, the Fermi data favor dark matter with many decay steps and a large DM mass while 
ATIC requires either a DM model with few decay steps or a significant local substructure in DM 
density distribution. 

Let us now turn to pulsars. In the calculation of electron and positron fiuxes, pulsars can be 
considered as point-like instantaneous sources [28], Q ~ 5{x — XQ)5{t — to). The main reason is 
that the typical propagation time 100 kyr) is much larger than the characteristic time scale 
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FIG. 2: The concordance model corresponding to the triangle in Figure [T] with the parameters given in 
Table |l} The upper (lower) edge of the band is the best fit for ATIC (Fermi) + PAMELA points with 
energies > 10 GeV. HESS 2008 and HESS 2009 points fit the concordance model, if multiphed by a 
factor ~ 0.9 which may be related to contamination by diffuse gamma rays. 



when a pulsar loses most of its rotational energy and the electrons and positrons are released to 
the interstellar medium (ISM). 

For every pulsar, we will consider two energy scales [131 EH]: the cutoff in the injection 
spectrum of electrons and positrons from the pulsar into the ISM, -Emj.cut, which can be between 
few hundred GeVs and tens of TeVs (see, e.g., [29] and references in [28]), and the cooling break 
i?br = obtained by integrating the energy losses E = —BqE^ during the propagation time 
t, where t is approximately the age of the pulsar. At E ^ 100 GeV the energy losses can be 



6 



estimated as 60 ~ 1.6GeV"^s-i ^ 5TeV"^Myr~^ [30]. At energies E > 100 GeV the coefficient 60 
slightly decreases with the energy [30j, since the Thompson approximation to inverse Compton 
scattering between electrons and the starlight becomes inapplicable. 

For pulsars with age t ^ 10 Myr, the cooling break is i^br ^ 20 GeV. In the ATNF catalog 
there are several hundred pulsars within 3 kpc from Earth and an age t < 10 Myr [31]. Below 
approximately 300 GeV, the corresponding flux is well approximated by the flux from a continuous 
distribution of pulsars in the Galactic plane [28]. An index 7^ 2 requires an index in the electron 
injection spectrum from pulsars 7inj ^ 2 [281 1221 [331 EH ESI ESI l37] . 

At energies E ^ 300 GeV, the flux receives contributions only from young pulsars t ^ 1 Myr 
within a smaller distance d ^ 1 kpc (since the propagation time is smaller). One may expect 
only of order ten such pulsars [31]. There is also a lower bound on the age of pulsars, a few tens 
of kyr, due to the fact that the electrons may still be trapped by the Pulsar Wind Nebulae. 

Depending on the relative value of -Einj.cut and -Ecut, there are two possibilities for the e~^e~ 
flux from pulsars: 

• -2'inj.cut ~ -^'cut; then the observed spectrum is naturally flat if we assume that the injection 
spectrum from pulsars is flat. This possibility is favored by the Fermi data [T3] . 

• -E'inj.cut ^ Ecut, then the cutoff in the observed spectrum is due to the coohng break which 
is much sharper than an exponential cutoff. One should also expect a series of steps due to 
consecutive cooling break cutoffs from different pulsars [2B]- This possibility is consistent 
with ATIC but may be in tension with Fermi. 

To summarize, both Fermi and ATIC require an additional source of electrons and positrons 
with an index 7^ 2 at low energies and a cutoff ii^cut ~ 500 GeV. However the presence of a 
bump at high energies in ATIC data contradicts the smooth spectrum of Fermi and HESS. The 
sources that produce featureless spectrum include the DM models with Mdm > 1 TeV and several 
steps in DM annihilation process as well as pulsars with an injection cutoff -Emj.cut ~ -E'cut- The 
sources that can produce ATIC bump include DM models with Mdm ~ -^cut and few decay steps 
or DM models with M^m ^ -E'cut and a signiflcant contribution from a local clump. Pulsars with 
-^'inj.cut ^ Ecnt may give a sharp cutoff in the observed e~^e~ flux and, possibly, step-like features 
due to a series of cooling break cutoffs from the youngest nearby pulsars. This possibility is in 
tension with the Fermi data but consistent with the ATIC. 
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